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(54) Title: COMMUNICATION SYSTEM NOISE CANCELLATION POWER SIGNAL CALCULATION TECHNIQUES 

m 

(57) Abstract: In order to enhance the quality of a communication signal derived from speech and noise, a filter (50) divides a 
2 communication signal into a plurality of frequency band signals. A calculator generates a plurality of power band signals each 

having a power band value and corresponding to one of the frequency band signals. The power band values are based on estimating, 
^ over a time period, the power of one of the frequency band signals. The time period is different for different ones of the frequency 
^ band signals. The power band values are used to calculate the weighting factors which are used to alter the frequency band signals 
^ that are combined to generate an improved conxmunication signal (170). 
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TITLE OF INVENTION 

COMMUNICATION SYSTEM NOISE CANCELLATION POWER SIGNAL 
CALCULATION TECHNIQUES 

BACKGROUND OF THE INVENTION 

This invention relates to communication system noise cancellation techniques, 
and more particularly relates to calculation of power signals used in such techniques. 

The need for speech quality enhancement in single-channel speech 
communication systems has increased in importance especially due to the tremendous 
growth in cellular telephony. Cellular telephones are operated often in the presence of 
high levels of environmental background noise, such as in moving vehicles. Such 
high levels of noise cause significant degradation of the speech quality at the far end 
receiver. In such circumstances, speech enhancement techniques may be employed to 
improve the quality of the received speech so as to increase customer satisfaction and 
encourage longer talk times. 

Most noise suppression systems utilize some variation of spectral subtraction. 
Figure LA shows an example of a typical prior noise suppression system that uses 
spectral subtraction. A spectral decomposition of the input nois\' speech-containing 
signal is first performed using the Filter Bank. The Filter Bank may be a bank of 
bandpass filters (such as in reference [1]. which is identified at the end of the 
description of the preferred embodiments). The Filter Bank decomposes the signal 
into separate frequency bands. For each band, power measurements are performed 
and continuously updated over time in the Noisy Signal Power & Noise Power 
Estimation block. These power measures are used to determine the signal-to-noise 
ratio (SNR) in each band. The Voice Activity Detector is used to distinguish periods 
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of speech activity from periods of silence. The noise power in each band is updated 
primarily during silence while the noisy signal power is tracked at all times. For each 
frequency band, a gain (attenuation) factor is computed based on the SNR of the band 
and is used to attenuate the signal in the band. Thus, each frequency band of the 
noisy input speech signal is attenuated based on its SNR. 

Figure IB illustrates another more sophisticated prior approach using an 
overall SNR level in addition to the individual SNR values to compute the gain 
factors for each band. (See also reference [2].) The overall SNR is estimated in the 
Overall SNR Estimation block. The gain factor computations for each band are 
performed in the Gain Computation block. The attenuation of the signals in different 
bands is accomplished by multiplying the signal in each band by the corresponding 
gain factor in the Gain Multiplication block. Low SNR bands are attenuated more 
than the high SNR bands. The amount of attenuation is also greater if the overall 
SNR is low. After the attenuation process, the signals in the different bands are 
recombined into a single, clean output signal. The resulting output signal will have an 
improved overall perceived quality. 

The decomposition of the input noisy speech-containing signal can also be 
performed using Fourier transform techniques or wavelet transform techniques. 
Figure 2 shows the use of discrete Fourier transform techniques (shown as the 
Windowing & FFT block). Here a block of input samples is transformed to the 
frequency domain. The magnitude of the complex frequency domain elements are 
attenuated based on the spectral subtraction principles described eariier. The phase of 
the complex frequency domain elements are left unchanged. The complex frequency 
domain elements are then transformed back to the time domain via an inverse discrete 
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Fourier transform in the EFFT block, producing the output signal. Instead of Fourier 
transform techniques, wavelet transform techniques may be used for decomposing the 
input signal. 

A Voice Activity Detector is part of many noise suppression systems. 
Generally, the power of the input signal is compared to a variable threshold level. 
Whenever the threshold is exceeded, speech is assumed to be present. Otherwise, the 
signal is assumed to contain only background noise. Such two-state voice activity 
detectors do not perform robustly under adverse conditions such as in cellular 
telephony environments. An example of a voice activity detector is described in 
reference [5]. 

\'arious implementations of noise suppression systems utilizing spectral 
subtraction differ mainly in the methods used for power estimation, gain factor 
determination, spectral decomposition of the input signal and voice activity detection. 
A broad overview of spectral subtraction techniques can be found in reference [3]. 
Several other approaches to speech enhancement, as well as spectral subtraction, are 
overviewed in reference [4]. 

Accurate noisy signal and noise power measures, which are performed for 
each frequency band, are critical to the performance of any adaptive noise 
cancellation system. In the past, inaccuracies in such power measures have limited 
the effectiveness of known noise cancellation systems. This invention addresses and 
provides one solution tor such problems, 

BRIEF SUMMARY OF THE INVENTION 
.A preferred embodiment of the invention is useful in a communication system 
for processing a communication signal derived from speech and noise. The preferred 
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embodiment can enhance the quality of the communication signal. In order to achieve 
this result, the communication signal is divided into a plurality of frequency band 
signals, preferably by a filter or by a digital signal processor. A plurality of power 
band signals each having a power band value and corresponding to one of the 
frequency band signals are generated. Each of the power band values is based on 
estimating over a time period the power of one of the frequency band signals, and the 
time period is different for at least two of the frequency band signals. Weighting 
factors are calculated based at least in part on the power band values, and the 
frequency band signals are altered in response to the weighting factors to generate 
weighted frequency band signals. The weighted frequency band signals are combined 
to generate a communication signal with enhanced quality. The foregoing signal 
generations and calculations preferably are accomplished with a calculator. 

By using the foregoing techniques, the power measurements needed to 
improve communication signal quality can be made with a degree of ease and 
accuracy unattained by the known prior techniques. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Figures 1 A and IB are schematic block diagrams of known noise cancellation 
systems. 

Figure 2 is a schematic block diagram of another form of a known noise 
cancellation system. 

Figure 3 is a functional and schematic block diagram illustrating a preferred 
form of adaptive noise cancellation system made in accordance with the invention. 
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Figure 4 is a schematic block diagram illustrating one embodiment of the 
invention implemented by a digital signal processor. 

Figure 5 is graph of relative noise ratio versus weight illustrating a preferred 
assignment of weight for various ranges of values of relative noise ratios. 
5 Figure 6 is a graph plotting power versus Hz illustrating a typical power 

spectral density of background noise recorded from a cellular telephone in a moving 
vehicle. 

Figure 7 is a curve plotting Hz versus weight obtained from a preferred form 
of adaptive weighting function in accordance with the invention. 
10 Figure 8 is a graph plotting Hz versus weight for a family of weighting curves 

calculated according to a preferred embodiment of the invention. 

Figure 9 is a graph plotting Hz versus decibels of the broad spectral shape of a 
typical voiced speech segment. 

Figure 10 is a graph plotting Hz versus decibels of the broad spectral shape of 
15 a typical unvoiced speech segment. 

Figure 1 1 is a graph plotting Hz versus decibels of perceptual spectral 
weighting curves for ko=25. 

Figure 12 is a graph plotting Hz versus decibels of perceptual spectral 
weighting curves for ko=38. 
20 Figure 13 is a graph plotting Hz versus decibels of perceptual spectral 

weighting curves for ko=50. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
The preferred form of ANC system shown in Figure 3 is robust under adverse 
conditions often present in cellular telephony and packet voice networks. Such 
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adverse conditions include signal dropouts and fast changing background noise 
conditions with wide dynamic ranges. The Figure 3 embodiment focuses on attaining 
high perceptual quality in the processed speech signal under a wide variety of such 
channel impairments. 

The performance limitation imposed by commonly used two-state voice activity 
detection functions is overcome in the preferred embodiment by using a probabilistic 
speech presence measure. This new measure of speech is called the Speech Presence 
Measure (SPM), and it provides multiple signal activity states and allows more 
accurate handling of the input signal during different states. The SPM is capable of 
detecting signal dropouts as well as new environments. Dropouts are temporary 
losses of the signal that occur commonly in cellular telephony and in voice over 
packet networks. New environment detection is the ability to detect the start of new 
calls as well as sudden changes in the background noise environment of an ongoing 
call. The SPM can be beneficial to any noise reduction function, including the 
preferred embodiment of this invention. 

Accurate noisy signal and noise power measures, which are performed for each 
frequency band, improve the performance of the preferred embodiment. The 
measurement for each band is optimized based on its frequency and the state 
information from the SPM. The frequency dependence is due to the optimization of 
power measurement lime constants based on the statistical distribution of power 
across the spectrum in typical speech and environmental background noise. 
Furthermore, this spectrally based optimization of the power measures has taken into 
consideration the non-linear nature of the human auditory system. The SPM state 
information provides additional information for the optimization of the time constants 
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as well as ensuring stability and speed of the power measurements under adverse 
conditions. For instance, the indication of a new environment by the SPM allows the 
fast reaction of the power measures to the new environment. 

According to the preferred embodiment, significant enhancements to perceived 
quality, especially under severe noise conditions, are achieved via three novel spectral 
weighting functions. The weighting functions are based on (1) the overall noise-to- 
signal ratio (NSR), (2) the relative noise ratio, and (3) a perceptual spectral weighting 
model. The first function is based on the fact that over-suppression under heavier 
overall noise conditions provide better perceived quality. The second function utilizes 
the noise contribution of a band relative to the overall noise to appropriately weight 
the band, hence providing a fine structure to the spectral weighting. The third 
weighting function is based on a model of the power-frequency relationship in typical 
environmental background noise. The power and frequency are approximately 
inversely related, from which the name of the model is derived. The inverse spectral 
weighting model parameters can be adapted to match the actual environment of an 
ongoing call. The weights are conveniently applied to the NSR values computed for 
each frequency band; although, such weighting could be applied to other parameters 
with appropriate modifications just as well. Funhermore, since the weighting 
functions are independent, only some or all the functions can be jointly utilized. 

The prefeiTed embodiment preserves the natural spectral shape of the speech 
signal which is important to perceived speech quality. This is attained by careful 
spectrally interdependent gain adjustment achieved through the attenuation factors. 
An additional advantage of such spectrally interdependent gain adjustment is the 
variance reduction of the attenuation factors. 
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Referring to Figure 3, a preferred form of adaptive noise cancellation system 
10 made in accordance with the invention comprises an input voice channel 20 
transmitting a communication signal comprising a plurality of frequency bands 
derived from speech and noise to an input terminal 22. A speech signal component of 
the communication signal is due to speech and a noise signal component of the 
communication signal is due to noise. 

A filter function 50 filters the communication signal into a plurahty of 
frequency band signals on a signal path 51. A DTMF tone detection function 60 and 
a speech presence measure funcuon 70 also receive the communication signal on 
input channel 20. The frequency band signals on path 51 are processed by a noisy 
signal power and noise power estimation function 80 to produce various forms of 
power signals. 

The power signals provide inputs to an perceptual spectral weighting function 
90, a relative noise ratio based weighting function 100 and an overall noise to signal 
ratio based weighting function 1 10. Functions 90, 100 and 1 10 also receive inputs 
from speech presence measure function 70 which is an improved voice activity 
detector. Functions 90, 100 and 1 10 generate preferred forms of weighting signals 
having weighting factors for each of the frequency bands generated by filter function 
50. The weighting signals provide inputs to a noise to signal ratio computation and 
weighting function 120 which multiplies the weighting factors from functions 90. 100 
and 110 for each frequency band together and computes an NSR value for each 
frequency band signal generated by the filter function 50. Some of the power signals 
calculated by function 80 also provide inputs to function 120 for calculating the NSR 
value. 
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Based on the combined weighting values and NSR value input from function 
120, a gain computation and interdependent gain adjustment function 130 calculates 
preferred forms of initial gain signals and preferred forms of modified gain signals 
with initial and modified gain values for each of the frequency bands and modifies the 
initial gain values for each frequency band by, for example, smoothing so as to reduce 
the variance of the gain. The value of the modified gain signal for each frequency 
band generated by function 130 is multiplied by the value of every sample of the 
frequency band signal in a gain multiplication function 140 to generate preferred 
forms of weighted frequency band signals. The weighted frequency band signals are 
summed in a combiner function 160 to generate a communication signal which is 
transmitted through an output terminal 172 to a channel 170 with enhanced quality. A 
DTMF tone extension or regeneration function 150 also can place a DTMF tone on 
channel 170 through the operation of combiner function 160. 

The function blocks shown in Figure 3 may be implemented by a variety of 
well known calculators, including one or more digital signal processors (DSP) 
including a program memory storing programs which are executed lo perform the 
functions associated with the blocks (described later in more detail) and a data 
memory for storing the variables and other data described in connection with the 
blocks. One such embodiment is shown in Figure 4 which illustrates a calculator in 
the form of a digital signal processor 12 which communicates with a memory 14 over 
a bus 16. Processor 12 performs each of the functions identified in connection with 
the blocks of Figure 3. Altemativeiy, any of the function blocks may be implemented 
by dedicated hardware implemented by application specific integrated circuits 
(ASICs), including memory, which are well known in the art. Of course, a 
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combination of one or more DSPs and one or more ASICs also may be used to 
implement the preferred embodiment. Thus, Figure 3 also illustrates an ANC 10 
comprising a separate ASIC for each block capable of performing the function 
indicated by the block. 
Filtering 

In typical telephony applications, the noisy speech-containing input signal on 
channel 20 occupies a 4kH2 bandwidth. This communication signal may be 
spectrally decomposed by filter 50 using a filter bank or other means for dividing the 
communication signal into a plurality of frequency band signals. For example, the 
filter function could be implemented with block-processing methods, such as a Fast 
Fourier Transform (FFT). In the case of an FFT implementation of filter function 50. 
the resulting frequency band signals typically represent a magnitude value (or its 
square) and a phase value. The techniques disclosed in this specification typically are 
applied to the magnitude values of the frequency band signals. Filter 50 decomposes 
the input signal into :V frequency band signals representing N frequency bands on 
path 51. The input to filter 50 will be denoted x{u) while the output of the k'" filter 
in the filter 50 will be denoted .v., (n) , where u is the sample time. 

The input. .v(;0 . to filter 50 is high-pass filtered to remove DC components by 
conventional means not shown. 

Gain Computation 

We first will discuss one form of gain computation. Later, we will discuss an 
interdependent gain adjustment technique. The gain (or attenuation) factor for the k''" 
frequency band is computed by function 130 once every T samples as 
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l'W,{n)NSR,in) , /i = 0,r,2r,... 

G,(/i-i) , 72 = 1,2 r-i,r + L,..,2r-L... 

A suitable value for T is 10 when the sampling rate is 8kHz. The gain factor will 
range between a small positive value, a , and 1 because the weighted NSR values are 
limited to lie in the range [0,1- £ ]. Setting the lower limit of the gain to e reduces the 
effects of ''musical noise'' (described in reference [2]) and permits limited background 
signal transparency. In the preferred embodiment, 8 is set to 0,05. The weighting 
factor, W^, (n ) , is used for over-suppression and under-suppression purposes of the 

signal in the k'^' frequency band. The overall weighting factor is computed by 
function 1 20 as 

W,(/7) = H,(72>v(«)w,(n) (2) 

where in) is the weight factor or value based on overall NSR as calculated by 

function 110, w^, (/?) is the weight factor or value based on the relative noise ratio 

weighting as calculated by function 100, and v. (/i)is the weight factor or value based 

on perceptual spectral weighting as calculated by function 90. As previously 
described, each of the weight factors may be used separately or in various 
combinations. 

Gain Multiplication 

The attenuation of the signal \\(n) from the k'^ frequency band is achieved 
by function 140 by multiplying \\(nj by its corresponding gain factor, G,.[u) , every 
sample to generate weighted frequency band signals. Combiner 160 sums the 
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resulting attenuated signals, y(n) , to generate the enhanced output signal on channel 
170. This can be expressed mathematically as: 

y{n) = ^G,in)x,{n) (3) 

Power Estimation 

The operations of noisy signal power and noise power estimation function 80 
include the calculation of power estimates and generating preferred forms of 
corresponding power band signals having power band values as identified in Table 1 
below. The power, P(n) sample 72, of a discrete-time signal u(n)y is estimated 
approximately by either (a) lowpass filtering the full-wave rectified signal or (b) 
lowpass filtering an even power of the signal such as the square of the signal. A first 
order IIR filter can be usedlfor the lowpass filter for both cases as follows: 

Pf/i) = l3P(n-l)-^a\u(n)\ f4a) 

P(n) = (3P{fi-l) + a[iiin)y- (4b) 

The lowpass filtering of the full-wave rectified signal or an even power of a signal is 
an averaging process. The power estimation (e.g., averaging) has an effective time 
window or time period during which the filter coefficients are large, whereas outside 
this window, the coefficients are close to zero. The coefficients of the lowpass filter 
determine the size of this window or lime period. Thus, the power estimation (e.g.. 
averaging) over different effective window sizes or time periods can be achieved by 
using different filter coefficients. When the rate of averaging is said to be increased, 
it is meant that a shoner time period is used. By using a shoner time period, the 
power estimates react more quickly to the newer samples, and '"forget" the effect of 
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older samples more readily. When the rate of averaging is said to be reduced, it is 
meant that a longer time period is used. 



The first order IIR filter has the following transfer function 

l-jSz 



The DC gain of this filter is //(I) = ^he coefficient, )8 , is a decay constant. 

The decay constant represents how long it would take for the present (non-zero) value 
of the power to decay to a small fraction of the present value if the input is zero, i.e. 
u{n) = 0 . If the decay constant, ^ . is close to unity, then it will take a longer time 
for the power value to decay. If /3 is close to zero, then it will take a shorter time for 
the power value to decay. Thus, the decay constant also represents how fast the old 
power value is forgotten and how quickly the power of the newer input samples is 
incorporated. Thus, larger values of jS result in longer effective averaging windows 
or time periods. 

Depending on the signal of interest, effectively averaging over a shorter or 
longer time period may be appropriate for power estimation. Speech power, which 
has a rapidly changing profile, would be suitably estimated using a smaller . Noise 
can be considered stationary for longer periods of time than speech. Noise power 
would be more accurately estimated by using a longer averaging window (large ^ ). 

The preferred form of power estimation significantly reduces computational 
complexity by undersampling the input signal for power estimation purposes. This 
means that only one sample out of every T samples is used for updating the power 
P{n) in (4). Between these updates, the power estimate is held constant. This 
procedure can be mathematically expressed as 
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_ J /8P(« - 1) + a I uin) I , « = 0,27,37,... 
^"^ 1P(«-1) , ,i = l,2,...r-l.r + 1....2T-l.... 

Such first order lowpass IIR filters may be used for estimation of the various power 
measures listed in the Table 1 below: 



Variable 


Description 




Overall noisy signal power 




Overall background noise power 




Noisy signal power in the k'^ frequency 
band. 




Noise power in the k'^ frequency band. 




Short-term overall noisy signal power in 
the first formant 




Long-term overall noisy signal power in 
the first formant 



Function 80 generates a signal for each of the foregoing Variables. Each of the 
signals in Table 1 is calculated using the estimations described in this Power 
Estimation section. The Speech Presence Measure, which will be discussed later, 
utilizes short-term and long-term power measures in the first formant region. To 
perform the first formant power measurements, the input signal. .v(n) , is lowpass 

filtered usina an IIR filter H{z} = ^^^i- "^^o-^ ^j^^ preferred 

implementation, the filter has a cut-off frequency at 850Hz and has coefficients 
bQ = 0.1027 . /?! = 0.2053, = -0.9754 and = 0.4103. Denoting the output of 
this filter as .v.. (// ) , the shon-term and long-term first formant power measures can 
be obtained as follows: 
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Pu,.LT (") = Pu..iTA Pu,,LT (" " D + 0!,,,,^, {n)\ 

= Pu.,LT.iPu,.LT (" - 1) + «,„.tr.:k/oK. («)| 



if ^u,..r(")<^u,.5r(«) 
and DROPOUT = 0 
if P,„.,^(;z)>/>,,^(«) 
and DROPOUT = 0 
= />,„ (n - 1) if DROPOUT = 1 

DROPOUT in (8) will be explained later. The time constants used in the above 
difference equations are the same as those described in (6) and are tabulated below: 



(8) 



Time Constant 


Value 




1/16000 




15999/16000 




1/256 




255/256 




1/128 


Pu.sr 


127/128 1 



One effect of these time constants is that the short term first formant power measure is 
effectively averaged over a shorter time period than the long term first formant power 
measure. These time constants are examples of the parameters used to analyze a 
communication signal and enhance its quality. 



10 



Noise-to-Signal Ratio (NSR) Estimation 

Regarding overall NSR based weighting function 110, the overall NSR, 
NSR^,^.^^^ii{n) at sample n . is defined as 

The overall NSR is used to influence the amount of over-suppression of the signal in 
each frequency band and will be discussed later. The NSR for the k'" frequency band 
may be computed as 



(9) 
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(n) 

MSR,{n)=^^^ (10) 

Those skilled in the an recognize that other algorithms may be used to compute the 
NSR values instead of expression (10). 
Speech Presence Measure (SPM) 

Speech presence measure (SPM) 70 may utilize any known DTMF detection 
method if DTMF tone extension or regeneration functions 150 are to be performed. 
In the preferred embodiment, the DTMF flag will be 1 when DTMF activity is 
detected and 0 otherwise. If DTMF tone extension or regeneration is unnecessary, 
then the following can be understood by always assuming that DTMF=0. 

SPM 70 pnmarily performs a measure of the likelihood that the signal activity is 
due to the presence of speech. This can be quantized to a discrete number of decision 
levels depending on the applicanon. In the preferred embodiment, we use five levels. 
The SPM performs its decision based on the DTMF flag and the LEVEL value. The 
DTMF flag has been described previously. The LEVEL value will be described 
shoniy. The decisions, as quantized, are tabulated below. The lower four decisions 
(Silence to High Speech) will be referred to as SPM decisions. 



Table 1: Joint Speech Presence Measure and DTMF Activity decisions 



DTMF 


LEVEL 


Decision 


1 


X 


DTMF Activity Present 


0 


0 


Silence Probability 


0 1 1 ! Low Speech Probability 


0 i 2 


Medium Speech Probability 


0 ! 3 


High Speech Probability 
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In addition to the above multi-level decisions, the SPM also outputs two flags or 
signals, DROPOUT and NEWENV, which will be described in the following 
sections. 

Power Measurement in the SPM 

The novel multi-level decisions made by the SPM are achieved by using a 
speech likelihood related comparison signal and multiple variable thresholds. In our 
preferred embodiment, we derive such a speech likelihood related comparison signal 
by comparing the values of the first formant short-term noisy signal power estimate, 
Pjsisii^h and the first formant long-term noisy signal power esdmate, Phuiinl 
Multiple comparisons are performed using expressions involving PLu.si^n) and 
P]st,Lj(^^) as given in the prefeired embodiment of equation (11) below. The result of 
these comparisons is used to update the speech likeHhood related comparison signal. 
In our preferred embodiment, the speech likelihood related comparison signal is a 
hangover counter, h^.^ . Each of the inequalities involving Pjsr.siijV and Pjsuii^) tises 
different scaling values (i.e. the fX. *s). They also possibly may use different additive 
constants, although we use Po=2 for all of them. 

The hangover counter, li^^ , can be assigned a variable hangover period that is 

updated every sample based on multiple threshold levels, which, in the preferred 
embodiment, have been limited to 3 levels as follows: 

= max [/i^^, - I] if /^.,,,,(«) > A^:^u./.7^"'*- ^. 
= max[/^,,,,,/^,, -1] if P,,,st(^')> .^A..LT^^^'^^P^ 
- max[0,/i,^ - 1] otherwise 

^'here h^^. > /?^^^_, > /;^^, and > > . 
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Suitable values for the maximum values of li^,^ are /z^^ =2000, h^^ . =1400 and 
^max,i = 800 . Suitable scaling values for the threshold comparison factors are 
^3 = 3.0, /i-, = 2.0 and jUj = L6. The choice of these scaling values are based on the 
desire to provide longer hangover periods following higher power speech segments. 
Thus, the inequalities of (11) determine whether Pistsii^) exceeds PistxM by more 
than a predetermined factor. Therefore, //^.^ represents a preferred form of 

comparison signal, resulting from the comparisons defined in (11) and having a value 
representing differing degrees of likelihood that a portion of the input communication 
signal results from at least some speech. 

Since longer hangover periods are assigned for higher power signal segments, 
the hangover period length can be considered as a measure that is directly 
proportional to the probability of speech presence. Since the SPM decision is required 
to reflect the likelihood that the signal activity is due to the presence of speech, and 
the SPM decision is based partly on the LEVEL value according to Table 1. we 
determine the value for LEVEL based on the hangover counter as tabulated below. 



Condition 


Decision 




LEVEL=3 1 




LEVEL=2 1 




LEVEL=l 


= 0 


LEVEL=0 



SPM 70 generates a preferred form of a speech likelihood signal having \ alues 
corresponding to LEVELs 0-3. Thus, LEVEL depends indirectly on the power 
measures and represents varying likelihood that the input communication signal 
results from at least some speech. Basing LEVEL on the hangover counter is 
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advantageous because a certain amount of hysterisis is provided. That is, once the 
count enters one of the ranges defined in the preceding table, the count is constrained 
to stay in the range for variable periods of lime. This hysterisis prevents the LEVEL 
value and hence the SPM decision from changing too often due to momentary 
changes in the signal power. If LEVEL were based solely on the power measures, the 
SPM decision would tend to flutter between adjacent levels when the power measures 
lie near decision boundaries. 

Dropout Detection in the SPM 

Another novel feature of the SPM is the ability to detect 'dropouts' in the 
signal. A dropout is a situation where the input signal power has a defined attribute, 
such as suddenly dropping to a very low level or even zero for short durations of time 
(usually less than a second). Such dropouts are often experienced especially in a 
cellular telephony environment. For example, dropouts can occur due to loss of 
speech frames in cellular telephony or due to the user moving from a noisy 
environment to a quiet environment suddenly. During dropouts, the ANC system 
operates differently as will be explained later. 

Dropout detection is incorporated into the SPM. Equation (8) shows the use of 
a DROPOUT signal in the long-term (noise) power measure. During dropouts, the 
adaptation of the long-term power for the SPM is stopped or slowed significantly. 
This prevents the long-term power measure from being reduced drastically during 
dropouts, which could potentially lead to incorrect speech presence measures later. 

The SPM dropout detection utilizes the DROPOUT signal or flag and a 
counter, c . The counter is updated as follows every sample time. 
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Condition 


Decision/Action 




=0 

^ dropout 




Increment 



The following table shows how DROPOUT should be updated. 



Condition 


Decision/Action 




DROPOlJr=l 


Otherwise 


DROPOUT=0 



As shown in the foregoing table, the attribute of c^^^^^^ determines at least in part the 
condition of the DROPOUT signal. A suitable value for the power threshold 
comparison factor, yU^^^p^,, , is 0.2. Suitable values for and c. are = 4000 and 
c, = 8000 . which correspond to 0.5 and 1 second, respectively. The logic presented 
here prevents the SPM from indicating the dropout condition for more than 
samples. 

Limiting of Long-term (Noise) Power Measure in the SPM 

In addition to the above enhancements to the long-term (noise) power 
measure, P^^^ {n) , it is further consu-ained from exceeding a cenain threshold, 

Put LT max ' i-^' ^^'^^ ^i.v/.L7 computcd according lo equation (7) is greater 
than Pi,,.^r.max ^ ^hen we set P,,,^,rW = Put,LT,n^ • This enhancement to the long-term 
power measure makes the SPM more robust as it will not be able to rise to the level of 
the shon-term power measure in the case of a long and continuous period of loud 
speech. This prevents the SPM from providing an incorrect speech presence measure 
in such situations. A suitable value for P,,,.^^.^, = 500/8159 assuming that the 
maximum absolute value of the input signal .v(;z) is normalized to unity. 



wo 01/73760 



PCT/USOl/06889 



New Environment Detection in the SPM 

At the beginning of a call, the background noise environment would not be 
known by ANC system 10. The background noise environment can also change 
suddenly when the user moves from a noisy environment to a quieter environment 
e.g. moving from a busy street to an indoor environment with windows and doors 
closed. In both these cases, it would be advantageous to adapt the noise power 
measures quickly for a short period of time. In order to indicate such changes in the 
environment, the SPM outputs a signal or flag called NEWENV to the ANC system. 

The detection of a new environment at the beginning of a call will depend on 
the system under question. Usually, there is some form of indication that a new call 
has been initiated. For instance, when there is no call on a particular line in some 
networks, an idle code may be transmitted. In such systems, a new call can be 
detected by checking for the absence of idle codes. Thus, the method for inferring 
that a new call has begun will depend on the particular system. 

In the preferred embodiment of the SPM. we use the flag NEWENV together 
with a counter c,,^.^^,,^. and a flag, OLDDROPOUT. The OLDDROPOUT flag 

contains the value of the DROPOLT from the previous sample time. 

A pitch estimator is used to monitor whether voiced speech is present in the 
input signal. If voiced speech is present, the pitch period (i.e., the inverse of pitch 
frequency) would be relatively steady over a period of about 20ms. If only 
background noise is present, then the pitch period would change in a random manner 
If a cellular handset is moved from a quiet room to a noisy outdoor environment, the 
input signal would be suddenly much louder and may be incorrectly detected as 



wo 01/73760 



22 



PCT/USOl/06889 



speech. The pitch detector can be used to avoid such incoirect detection and to set the 
new environment signal so that the new noise environment can be quickly measured. 

To implement this function, any of the numerous known pitch period 
estimation devices may be used, such as device 74 shown in Fig. 3. In our preferred 
implementation, the following method is used. Denoting K(n-T) as the pitch period 
estimate from T samples ago, and K(n) as the current pitch period estimate, if I^Tfnj- 
K(n-40)\>3. and \K(n-40)'K(n'80)\>3. and \K(n'80)'K(n'I20)\>3, then the pitch 
period is not steady and it is unfikely that the input signal contains voiced speech. If 
these conditions are true and yet the SPM says that LEVEL>1 which normally implies 
that significant speech is present, then it can be inferred that a sudden increase in the 
background noise has occurred. 



The following table specifies a method of updating NEWENV and 



Condition 


Decision/Action j 


Beginning of a new call or 


NEWENV=:i . 


( (6lDDR0P0UT=1) and (DROPOUT=0) ) or 


c ==0 i 

fmvenv ' 


{\K(n)-K(n-40)\>3 and \K(n-40)-K(n-80)\>3 and 
\K(n-80)-K(n-120}\>3 and LEVEL>I) 


i 

1 

! 


Not the beginning of a new call or 


No action 1 


OLDDROPOUT=0 or 




DR0P0UT=1 




c <c andNEWENV=l 


Increment c„ \ 


C = C 
wwenv neH'fmMnax 


NEWENV=0 1 




c =0 ! 

'tewenv \ 



In the above method, the NEWENV flag is set to 1 for a period of time specified by 
^neuw.max ' ^^^^^ which it js clcarcd. The NEWENV flag is set to 1 in response to 
various events or attributes: 



f 1) at the beginning of a new call: 
(2) at the end of a dropout period: 
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(3) in response to an increase in background noise (for example, the pitch 
detector 74 may reveal that a new high amplitude signal is not due to speech, but 
rather due to noise.); or 

(4) in response to a sudden decrease in background noise to a lower level of 
sufficient amplitude to avoid being a drop out condition. 

A suitable value for the c,,^,„, ^, is 2000 which corresponds to 0.25 seconds. 

Operation of the ANC System 

Referring to Figure 3. the multi-level SPM decision and the flags DROPOUT 
and NEWENV are generated on path 72 by SPM 70. With these signals, the ANC 
system is able to perform noise cancellation more effectively under adverse 
conditions. Funhermore, as previously described, the power measurement function 
has been significantly enhanced compared to prior known systems. Additionally, the 
three independent weighting functions carried out by functions 90, 100 and 110 can 
be used to achieve over-suppression or under-suppression. Finally, gain computation 
and interdependent gain adjustment function 130 offers enhanced performance. 

Use of Dropout Signals 

When the flag DROPOUT=l. the SPM 70 is indicating that there is a 
temporary loss of signal. Under such conditions, continuing the adaptation of the 
signal and noise power measures could result in poor behavior of a noise suppression 
system. One solution is to slow down the power measurements by using very long 
time constants. In the preferred embodiment, we freeze the adaptation of both signal 
and noise power measures for the individual frequency bands, i.e. we.set 
P; (n) = P^ (;i-l) and P; (;n = 7^; (;i - 1) when DR0P0UT=1. Since DROPOUT 



wo 01/73760 



24 



PCTAJSOl/06889 



remains at 1 only for a short time (at most 0.5 sec in our implementation), an 
erroneous dropout detection may only affect ANC system 10 momentarily. The 
improvement in speech quality gained by our robust dropout detection outweighs the 
low risk of incorrect detection. 

Use of New Environment Signals 

When the flag NEWENV=1, SPM 70 is indicating that there is a new 
environment due to either a new call or that it is a post-dropout environment. If there 
is no speech activity, i.e. the SPM indicates that there is silence, then it would be 
advantageous for the ANC system to measure the noise spectrum quickly. This quick 
reaction allows a shorter adaptation time for the ANC system to a new noise 
environment. Under normal operation, the time constants, a ^ Pi- > ^sed for the 
noise power measurements would be as given in Table 2 below. When NEWENV=1, 
we force the time constants to correspond to those specified for the Silence state in 
Table 2. The larger 3 values result in a fast adaptation to the background noise power. 
SPM 70 will only hold the NEWENV at 1 for a short period of time. Thus, the ANC 
system will automatically revert to using the normal Table 2 values after this time. 
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Table 2: Power measurement time constants 



SPM 
Decision 


Frequency Range 


Time Constants 






or* 


Pj 


Silence Probability 
LEVEL=0 


<800Hzor>2500Hz 


r/60 


i-r/6000 


0.533 


1-7/240 


800Hz to 2500Hz 


r/80 


1-7-/8000 


0.533 


1-7/240 


Low Speech 
Probability 

T PVPT —1 


<800Hz or >2500Hz 


r/120 


i-r/12000 


0.533 


1-7/240 


800Hz to 2500Hz 


r/160 


i-r/16000 


0.64 


l-T 1 200 


Medium Speech 
Probability 
LEVEL=2 


<800Hzor>2500Hz 


Noise power values 
remain substantially 
constant. 


0.64 


1-7/200 


800Hz to 2500Hz 


0.853 


1-7/150 


High Speech 
Probability 
LEVEL=3 


<800Hz or>2500Hz 


0.853 


1-7/150 


800Hz to 2500Hz 


1 


1-7/128 



10 



15 



Frequencv-Dependent and Speech Presence Measure-Based Time Constants for 
Power Measurement 

The noise and signal power measurements for the different frequency bands are given 
by 

'^^"^ IP^'(n-l) . /i = 1.2,...7-l.r + 1....2r-U. 

' \p',{n'Y) . /7 = u,...7-ir + i....:r-u. 

In the preferred embodiment, the time constants ^ y , ^5 , 0:^ and are based on 

both the frequency band and the SPM decisions. The frequency dependence will be 
explained first, followed by the dependence on the SPM decisions. 

The use of different time constants for power measurements in different 
frequency bands offers advantages. The power in frequency bands in the middle of 
the 4kHz speech bandwidth naturally tend to ha\*e higher average power levels and 
variance during speech than other bands. To track the faster variations, it is useful to 
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have relatively faster time constants for the signal power measures in this region. 
Relatively slower signal power time constants are suitable for the low and high 
frequency regions. The reverse is true for the noise power time constants, i.e. faster 
time constants in the low and high frequencies and slower time constants in the 
middle frequencies. We have discovered that it would be better to track at a higher 
speed the noise in regions where speech power is usually low. This results in an 
earlier suppression of noise especially at the end of speech bursts. 

In addition to the variation of time constants with frequency, the time 
constants are also based on the multi-level decisions of the SPM. In our preferred 
implementation of the SPM. there are four possible SPM decisions (i.e.. Silence, Low 
Speech, Medium Speech, High Speech). When the SPM decision is Silence, it would 
be beneficial to speed up the tracking of the noise in all the bands. When the SPM 
decision is Low Speech, the likelihood of speech is higher and the noise power 
measurements are slowed down accordingly. The likelihood of speech is considered 
too high in the remaining speech states and thus the noise power measurements are 
turned off in these states. In contrast to the noise power measurement, the time 
constants for the signal power measurements are modified so as to slow down the 
tracking when the likelihood of speech is low. This reduces the variance of the signal 
power measures during low speech levels and silent periods. This is especially 
beneficial during silent periods as it prevents shon-duration noise spikes from causing 
the gain factors to rise. 

In the preferred embodiment, we have selected the time constants as shown in 
Table 2 above. The DC gains of the IIR filters used for power measurements remain 
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fixed across all frequencies for simplicity in our preferred embodiment although this 
could be varied as well. 

Weighting based on Overall NSR 

In reference [2], it is explained that the perceived quality of speech is 
improved by over-suppression of frequency bands based on the overall SNR. In the 
preferred embodiment, over-suppression is achieved by weighting the NSR according 
to (2) using the weight, iii. (n) , given by 

Here, we have limited the weight to range from 0.5 to 1.5. This weight computation 
may be performed slower than the sampling rate for economical reasons. A suitable 
update rate is once per 2T samples. 

Weighting Based on Relative Noise Ratios 

We have discovered that improved noise cancellation results from vveightmg 
based on relative noise ratios. According to the preferred embodiment, the weighting, 
denoted by uv , based on the values of noise power signals in each frequency band, 
has a nominal value of unity for all frequency bands. This weight will be higher for a 
frequency band that contributes relatively more to the total noise than other bands. 
Thus, greater suppression is achieved in bands that have relatively more noise. For 
bands chat contribute little to the overall noise, the weight is reduced below unitv to 
reduce the amount of suppression. This is especially imponant when both the speech 
and noise power m a band are very low and of the same order. In the past, in such 
Siiuatiuns.. pouer has been scv^-reiy suppressed, u hich has resulted in noilow 
sounding speech. However, with this weighting function, the amount of suppression 
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is reduced, preserving the richness of the signal, especially in the high frequency 
region. 

There are many ways to determine suitable values for vv^. , First, we note that 
the average background noise power is the sum of the background noise powers in N 
frequency bands divided by the N frequency bands and is represented by Ps^in)! N , 
The relative noise ratio in a frequency band can be defined as 

The goal is to assign a higher weight for a band when the ratio, R,. (n) . for thai 

band is high, and lower weights when the ratio is low. In the preferred embodiment, 
we assign these weights as shown in Figure 5, where the weights are allowed to range 
between 0.5 and 2. To save on computational time and cost, we perform the update of 
(15) once per 2T samples. Function 80 (Figure 3) generates preferred forms of band 
power signals corresponding to the terms on the right side of equation (15) and 
function 100 generates preferred forms of weighting signals with weighting values 
corresponding to the term on the left side of equation (15). 

If an approximate knowledge of the nature of the environmental noise is 
known, then the RNR weighting technique can be extended to incorporate this 
knowledge. Figure 6 shows the typical power spectral density of background noise 
recorded from a cellular telephone in a moving vehicle. Typical environmental 
background noise has a power spectrum that corresponds to pink or brown noise. 
(Pink noise has power inversely proportional to the frequency. Brown noise has 
poAcr ;n\erseiy proponion.il ro ihe square of ihe :'requer;c\. ; 3;.:sed on this 
approximate knowledge of the relative noise ratio profile across the frequency bands. 
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the perceived quality of speech is improved by weighting the lower frequencies more 
heavily so that greater suppression is achieved at these frequencies. 

We take advantage of the knowledge of the typical noise power spectrum 
profile (or equivalently, the RNR profile) to obtain an adaptive weighting function. In 
general, the weight, w- for a particular frequency, / , can be modeled as a function 

of frequency in many ways. One such model is 

=b{f~f,)'^c (16) 
This model has three parameters [b.f^.c]. An example of a weighting curve 

obtained from this model is shown in Figure 1 for b - 5.6x10"** , = 3000 and 
c = 0.5 , 

The Figure 7 curve varies monotonically with decreasing values of weight from 0 Hz 
to about 3000 Hz, and also vanes monotonically with increasing values of weight 
from about 3000 Hz to about 4000 Hz. In practice, we could use the frequency band 
index, k , corresponding to the actual frequency / . This provides the following 
practical and efficient model with parameters { b.k^x ): 

=b(k~k^)- (17) 
In general, the ideal weights, w^. , may be obtained as a function of the measured noise 

power estimates, , at each frequency band as follows: 

=min(L ^— ) (IS) 

max{P;} 

jt 

Basically, the ideal weights are equal to the noise power measures normalized by the 
largest noise power measure. In general, the normalized power of a noise component 
in a panicular frequency band is defined as a ratio of the power of the noise 
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component in that frequency band and a function of some or all of the powers of the 
noise components in the frequency band or outside the frequency band. Equations 
(15) and (18) are examples of such normalized power of a noise component. In case 
all the power values are zero, the ideal weight is set to unity. This ideal weight is 
actually an alternative definition of RNR. We have discovered that noise cancellation 
can be improved by providing weighting which at least approximates normalized 
power of the noise signal component of the input communication signal. In the 
preferred embodiment, the normalized power may be calculated according to (18). 
Accordingly, function 100 CFigure 3) may generate a preferred form of weighting 
signals having weighting values approximating equation (18). 

The approximate model in (17) attempts to mimic the ideal weights computed 
using (18). To obtain the model parameters {b^k^^c least-squares approach may 
be used. An efficient way to perform this is to use the method of steepest descent to 
adapt the model parameters { ,c }. 

We derive here the general method of adapting the model parameters using the 
steepest descent technique. First, the total squared error between the weights 
generated by the model and the ideal weights is defined for each frequency band as 
follows: 



ail ;t 

Taking the panial derivative of the total squared error, e' , with respect to each of the 
model parameters in turn and droppmg constant terms, we obtain 




(19) 
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-^ = -^[bik-k,y-+c- w, }?(k -k,) (21) 
= Ykk-k,y-+c-w,] (22) 

all it 

Denoting the model parameters and the error at the n'^ sample time as { b^.k^^ , c„ } 
and (k) , respectively, the model parameters at the (n -h I)''' sample can be estimated 



^0 ^li 

del 
dc 



as 



^n.l =^n-A,-^ (23) 

J J . de' 

=^0.;, — (24) 

dc 



Here { A^.A,. .A. } are appropriate step-size parameters. The model definition in (17) 
can then be used to obtain the weights for use in noise suppression, as well as being 
used for the next iteration of the algorithm. The iterations may be performed everv' 
sample time or slower, if desired, for economy. 

We have described the alternative preferred RNR weight adaptation technique 
above. The weights obtained by this technique can be used to directly multiply the 
corresponding NSR values. These are then used to compute the gain factors for 
attenuation of the respective frequency bands. 

In another embodiment, the weights are adapted efficiently using a simpler 
adaptation technique for economical reasons. We fi.x the value of the weishtins 
modei parameter lo a:,, = 36 which corresponds to = 2SS0Hz in t i6). 
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Furthermore, we set the model parameter at sample time n to be a function of fc^j 
and the remaining model parameter as follows: 

^„=^ (26) 
Equation (26) is obtained by setting k=0 and u;. = 1 in (17). We adapt only c„ to 
determine the curvature of the relative noise ratio weighung curve. The range of is 
restricted to [0.1.1.0]. Several weighting curves corresponding to these specifications 
are shown in Figure 8. Lower values of c^, correspond to the lower curves. When 
= 1 , no spectral weighting is performed as shown in the uppermost line. For all 
other values of , the curves vary monotonically in the same manner described in 
connection with Figure 7. The greatest amount of curvature is obtained when 
c. =0.1 as shown in the lowest curve. The applicants have found it advantageous to 

arrange the weighting \ alues so that they var\' monotonically between two 
frequencies separated by a factor of 2 (e.g., the weighting values vary monotonically 
between 1000-2000 Hz and/or between 1500-3000 Hz), 

The deteiTnination of c^, is performed by comparing the total noise power in 

the lower half of the signal bandwidth to the total noise power in the upper half. We 
define the total noise power in the lower and upper half bands as: 

^.../.;,n..(''»= S^v('?) (27) 

Alternatively, lowpass and highpass filter could be used to filter xin) followed bv 
appropriate power measurement using (6) to obtain these noise powers. In our filter 
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bank implementation, /:e {3,4,...,42} and hence f;^^^,^ = {3,4,...22} and 
^upper - {23,24,. .,42} . Although these power measures may be updated every sample, 
they are updated once every 2T samples for econoniical reasons. Hence the value of 
needs to be updated only as often as the power measures. It is defined as follows: 



c,, = max 



nun 



1.0 



,0.1 



^UHol .lower ) 

The min and max functions restrict c^, to lie within [O.Ll.O]. 

According to another embodiment, a curve, such as Figure 7, could be stored 
as a weighting signal or table in memory 14 and used as static weighting values for 
each of the frequency band signals generated by filter 50. The curve could vary 
monotonically, as previously explained, or could vary according to the estimated 
spectral shape of noise or the estimated overall noise power. P^^, (n) ,as explained in 
the next paragraphs. 

Alternatively, the power spectral density shown in Figure 6 could be thought 
ot as defining the spectral shape of the noise component of the communication signal 
received on channel 20. The value of c is altered according to the spectral shape in 
order to determine the value of in equation (17). Spectral shape depends on the 
power of the noise component of the communication signal received on channel 20. 
As shown in equations (12) and (13). power is measured using time constants a{. and 
^ v which vary according to the likelihood of speech as shown in Table 2. Thus, the 
weighting values determined according to the spectral shape of the noise component 
ot the communication signal on channel 20 are derived in pan from the iii<eiihood tiiac 
the communication signal is derived at least in pan from speech. 



(29) 
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According to another embodiment, the weighting values could be determined 
from the overall background noise power. In this embodiment, the value of c in 
equation (17) is determined by the value of P^^ (n) . 

In general, according to the preceding paragraphs, the weighting values may 
vary in accordance with at least an approximation of one or more characteristics (e.g., 
spectral shape of noise or overall background power) of the noise signal component of 
the communication signal on channel 20. 
Perceptual Spectral Weighting 

We have discovered that improved noise cancellation results from perceptual 
spectral weighting (PSW) in which different frequency bands are weighted differently 
based on their perceptual importance. Heavier weighting results in greater 
suppression in a frequency band. For a given SNR (or NSR), frequency bands where 
speech signals are more important to the perceptual quality are weighted less and 
hence suppressed less. Without such weighting, noisy speech may sometimes sound 
'hollow' after noise reduction. Hollow sound has been a problem in previous noise 
reduction techniques because these systems had a tendency to oversuppress the 
perceptually important parts of speech. Such oversuppression was partly due to not 
taking into account the perceptually important spectral interdependence of the speech 
signal. 

The perceptual importance of different frequency bands change depending on 
characteristics of the frequency distribution of the speech component of the 
communication signal being processed. Determining perceptual importance from 
such characteristics may be accomplished by a variety of methods. For example, the 
characteristics may be determined by the likelihood that a communication signal is 
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derived from speech. As explained previously, this type of classification can be 
implemented by using a speiech likelihood related signal, such as h^..^ , Assuming a 

signal was derived from speech, the type of signal can be further classified by 
determining whether the speech is voiced or unvoiced. Voiced speech results from 
vibration of vocal cords and is illustrated by utterance of a vowel sound. Unvoiced 
speech does not require vibration of vocal cords and is illustrated by utterance of a 
consonant sound. 

The broad spectral shapes of typical voiced and unvoiced speech segments are 
shown in Figures 9 and 10, respectively. Typically, the lOOOHz to 3000H2 regions 
contain most of the power in voiced speech. For unvoiced speech, the higher 
frequencies (>2500Hz) tend to have greater overall power than the lower frequencies. 
The weighting in the PSW technique is adapted to maximize the perceived quality as 
the speech spectrum changes. 

As in RNR weighting technique, the actual implementation of the perceptual 
spectral weighting may be performed directly on the gain factors for the individual 
frequency bands. .Another aliemaiive is to weight the power measures appropnately. 
In our preferred method, the weighting is incorporated into the NSR measures. 

The PSW technique may be implemented independently or in any combination 
with the overall NSR based weighting and RNR based weighting methods. In our 
preferred implementation, we implement PSW together with the other two techniques 
as given in equation i2). 

The weigiits in the PSW technique are selected lo vary beiween zei'o ana one. 
Larger weights correspond to greater suppression. The basic idea of PSW is to adapt 
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the weighting curve in response to changes in the characteristics of the frequency 
distribution of at least some components of the communication signal on channel 20. 
For example, the weighting curve may be changed as the speech spectrum changes 
when the speech signal transitions from one type of communication signal to another, 
e.g., from voiced to unvoiced and vice versa. In some embodiments, the wei^htin^^ 
curve may be adapted to changes in the speech component of the communication 
signal. The regions that are most critical to perceived quality (and which are usually 
oversuppressed when usmg previous methods) are weighted less so that they are 
suppressed less. However, if these perceptually important regions contain a 
significant amount of noise, then their weights will be adapted closer to one. 

Many weightmg models can be devised to achieve the PSW, In a manner similar 
to the RNR technique's weighting scheme given by equation (17), we utilize the 
practical and efficient model with parameters {b.k^.c}: 

V, =b(k-k^r +c (30) 
Here \\ is the weight tor frequency band L In this method, we will vary only 

and c. This weighting curve is generally U-shaped and has a minimum value of c at 

frequency band kf^. For simplicity, we fix the weight at )t=0 to unity. This gives the 

following equation for /; as a function of k^^ and c: 

h='-^ (31. 
The lowest weight frequency band, k^, is adapted based on the likelihood of 
speech being voiced or i:n\*o:ced. in our prefen'ed method. A is allowed lo be in i.he 
range [25.50], which corresponds to the frequency range [2000Hz. 4000H2]. Dunns 
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Strong voiced speech, it is desirable to have the U-shaped weighting curve to have 
the lov^est weight frequency band kQ to be near 2000Hz. This ensures that the 
midband frequencies are weighted less in general. During unvoiced speech, the 
lowest weight frequency band is placed closer to 4000Hz so that the mid to high 

frequencies are weighted less, since these frequencies contain most of the perceptually 
important parts of unvoiced speech. To achieve this, the lowest weight frequency 
band is varied with the speech likelihood related comparison signal which is the 

hangover counter, /i^.^ , in our preferred method. Recall that h^.,^ is always in the 

range [0,/i^^^ ;=2000]. Larger values of /j^.^ indicate higher likelihoods of speech and 

also indicate a higher likelihood of voiced speech. Thus, in our preferred method, the 
lowest weight frequency band is varied with the speech likelihood related comparison 
signal as follows: 

^•o ==L50-/^J80J . (32) 

Since k^ is an integer, the floor function [_.J is used for rounding. 
Next, the method for adapting the minimum weight c is presented. In one 

approach, the minimum weight c could be fixed to a small value such as 0.25. 

However, this would always keep the weights in the neighborhood of the lowest 

weight frequency band k^ at this minimum value even if there is a strong noise 

component in that neighborhood. This could possibly result in insufficient noise 
attenuation. Hence we use the novel concept of a regional NSR to adapt the minimum 
weight. 
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The regional NSR, NSR^^^.^^^^ {k),is defined with respect to the minimum weight 
frequency band and is given by: 

Basically, the regional NSR is the ratio of the noise power to the noisy signal 
power in a neighborhood of the minimum weight frequency band . In our preferred 
method, we use up to 5 bands centered at k^ as given in the above equation. 

In our preferred implementation, when the regional NSR is -15dB or lower, we set 
the minimum weight c to 0.25 (which is about 12dB). As the regional NSR 
approaches its maximum value of OdB. the minimum weight is increased towards 
unity. This can be achieved by adapting the minimum weight c at sample time n as 

^^i 0,25 . .m.,^,,;,,(>0< 0.1778 = -15^5 

|0.912.V5/?_,(//)-f 0.088 . 0.1778 <//5/?_,in)< I ^^^^ 

The V, curves are plotted for a range of values of c and k^ in Figures 1 1-13 to 
illustrate the flexibility that this technique provides in adapting the weighting curves. 
Regardless of A',,, the curves are flat when c=L which corresponds to the situation 
where the regional NSR is unity (OdB). The curves shown in Figures 11-13 have the 
same monotonic properties and may be stored in memory 14 as a weighting signal or 
table in the same manner previously described in connection with Figure 7. 

As can be seen from equation (32). processor 12 generates a control signal from 
the speech likelihood signal h^..^^ which represents a characteristic of the speech and 
noise components of the communication signal on channel 20. As previously 
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explained, the likelihood signal can also be used as a measure of whether the speech is 
voiced or unvoiced. Determining whether the speech is voiced or unvoiced can be 
accomplished by means other than the likelihood signal. Such means are known to 
those skilled in the field of communications. 

The characteristics of the frequency distribution of the speech component of the 
channel 20 signal needed for PSW also can be determined from the output of pitch 
estimator 74. In this embodiment, the pitch estimate is used as a control signal which 
indicates the characteristics of the frequency distribution of the speech component of 
the channel 20 signal needed for PSW. The pitch estimate, or to be more specific, the 
rate of change of the pitch, can be used to solve for ko in equation (32). A slow rate 
of change would correspond to smaller ko values, and vice versa. 

In one embodiment of PSW. the calculated weights for the different bands are 
based on an approximation of the broad spectral shape or envelope of the speech 
component of the communication signal on channel 20. More specifically, the 
calculated weighting curve has a generally inverse relationship to the broad spectral 
shape of the speech component of the channel 20 signal. An example of such an 
inverse relationship is to calculate the weighting curve to be inversely proponional to 
the speech spectrum, such that when the broad spectral shape of the speech spectrum 
is multiplied by the weighting curve, the resulting broad spectral shape is 
approximately flat or constant at all frequencies in the frequency bands of interest. 
This is different from the standard spectral subtraction weighting which is based on 
the noise-to-signal ratio of individual bands. In this embodiment of PSW. w e are 
taking into consiucraiion ihe enure speech signal (or a significant ponion or" in to 
determine the weighting curve for all the frequency bands. In spectral subtraction, the 



wo 01/73760 



40 



PCT/USOl/06889 



weights are determined based only on the individual bands. Even in a spectral 
subtraction implementation such as in Figure IB. only the overall SNR or NSR is 
considered but not the broad spectral shape. 
Computation of Broad Spectral Shape or Envelope of Speech 

There are many methods available to approximate the broad spectral shape of 
the speech component of the channel 20 signal. For instance, linear prediction 
analysis techniques, commonly used in speech coding, can be used to determine the 
spectral shape. 

Alternatively, if the noise and signal powers of individual frequency bands are 
tracked using equations such as (12) and (13), the speech spectrum power at the k"^ 
band can be estimated as j P,^(/?)- P' (n)l. Since the soal is to obtain the broad 

spectral shape, the total power, P/ (;?) , may be used to approximate the speech power 
in the band. This is reasonable since, when speech is present, the signal spectrum 
shape is usually dominated by the speech spectrum shape. The set of band power 
values together provide the broad spectral shape estimate or envelope estimate. The 
number of band power values in the set will vary depending on the desired accuracy 
of the estimate. Smoothing of these band power values using moving average 
techniques is also beneficial to remove jaggedness in the envelope estimate. 
Computation of Perceptual Spectral Weighting Curve 

After the broad spectral shape is approximated, the perceptual weighting curve 
may be determined to be in\ ersely propoaionai to the broad spectral shape 

upproximaiion. For insuii-wj. w' :/m is used us :he bi-oac specira! shape estimate 
the k"^ band, then the weight for the k"^ band, v, , may be determined as 
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in) = \f// in) , where i// is a predetermined value. In this embodiment, a set of j 
speech power values, such as a set of {n) values, is used as a control signal 
indicating the characteristics of the frequency distribution of the speech component of 
the channel 20 signal needed for PSW. By using the foregoing spectral shape 
estimate and weighting curve, the variation of the power signals used for the estimate 
is reduced across the N frequency bands. For instance, the spectrum shape of the 
speech component of the channel 20 signal is made more nearly flat across the N 
frequency bands, and the variation in the spectrum shape is reduced. 

For economical reasons, we use a parametric technique in our preferred 
implementation which also has the advantage that the weighting curve is always 
smooth across frequencies. We use a parametric weighting curve, i.e. the weighting 
curve is formed based on a few parameters that are adapted based on the spectral 
shape. The number of parameters is less than the number of weighting factors. The 
parametric weighting function in our economical implementation is given by the 
equation (30), which is a quadratic cur\*e with three parameters. 

Use of Weighting Functions 

Although we have implemented weighting functions based on overall NSR 
( Uf. ), perceptual spectral weighting ( v,. ) and relative noise ratio weighting ( >v. ) 
jointly, a noise cancellation system will benefit from the implementation of only one 
or various combinations of the functions. 

In our preferred embodiment, we implement the weighting on the NSR values 
for the different tVequency bands. One could implement these weighting functions 
jusi as well, after appropriate modifications, directly on the gain factors. 
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Alternatively, one could apply the weights directly to the power measures prior to 
computation of the noise-to-signal values or the gain factors. A further possibility is 
to perform the different weighting functions on different variables appropriately in the 
ANC system. Thus, the novel weighting techniques described are not restricted to 
specific implementations. 

Spectral Smoothing and Gain Variance Reduction Across Frequency Bands 

In some noise cancellation applications, the bandpass filters of the filter bank 
used to separate the speech signal into different frequency band components have little 
overiap. Specifically, the magnitude frequency response of one filter does not 
significantly overlap the magnitude frequency response of any other filter in the filter 
bank. This is also usually true for discrete Fourier or fast Fourier transform based 
implementations. In such cases, we have discovered that improved noise cancellation 
can be achieved by interdependent gain adjustment. Such adjustment is affected by 
smoothing of the input signal spectrum and reduction in variance of gain factors 
across. the frequency bands according to the techniques described below. The 
splitting of the speech signal into different frequency bands and applying 
independently determineci gain factors on each band can sometimes destroy the 
natural spectral shape of the speech signal. Smoothing the gain factors across the 
bands can help to preserve the natural spectral shape of the speech signal. 
Furthermore, it also reduces the variance of the gain factors. 

This smoothing of the gain factors, G,.(/7) (equation (1)), can be performed by 
modifying each of the initial gain factors as a function of at least two of the initial 
gain factors. The initial gain factors preferably are generated in the form of signals 
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with initial gain values in function block 130 (Figure 3) according to equation (1). 
According to the preferred embodiment, the initial gain factors or values are modified 
using a weighted moving average. The gain factors corresponding to the low and 
high values of k must be handled slightly differently to prevent edge effects. The 
initial gain factoid are modified by recalculating equation (1) in function 130 to a 
preferred form of modified gain signals having modified gain values or factors. Then 
the modified gain factors are used for gain multiplication by equation (3) in function 
block 140 (Figure 3). 

More specifically, we compute the modified gains by first computing a set of 
initial gain values. Gl (n). We then perform a moving average weighting of these 
initial gain factors with neighboring gain values to obtain a new set of gain values, 
Gj,. (n) . The modified gain values derived from the initial gain values is given by 

C,(/z)= jM,G;(/n (35) 

The M are the moving average coefficients tabulated below for our preferred 
embodiment. 



Ranee of k 
*— 


Moving Average Weighting 
Coefficients. 


First coefficient to j 
be multiplied with | 


k = 3 


0.95. 0.04. 0.01 


GUu) 1 


k = 4 


0.02. 0.95, 0.02. 0.01 


G'An) 1 


5<k<40 


0.005. 0.02. 0.95. 0.02. 0.005 




A- = 41 


0.01.0.02. 0.95.0.02 




A- = 42 1 


0.01.0.04. 0.95 


G'^(n) 1 



We have discovered that improved noise cancellation is possible with 
coelTicients selected from the foilowing ranges or \ alues. One of ih^ coerncienis i:, \:: 
the range of 10 to 50 times the value of the sum of the other coefficients. For 
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example, the coefficient 0.95 is in the range of 10 to 50 times the value of the sum of 
the other coefficients shown in each line of the preceding table. More specifically, the 
coefficient 0.95 is in the range from .90 to .98. The coefficient 0.05 is in the range .02 
to .09. 

In another embodiment, we compute the gain factor for a particular frequency 
band as a function not only of the corresponding noisy signal and noise powers, but 
also as a function of the neighboring noisy signal and noise powers. Recall equation 



\l-WAn)NSR..{u) . = 0.7.27.... 
G.(/z) = ^ 

[ G,(/2-l) . ;2 = 1.2 7 - IT + 1,.... 27-1,... 

In this equation, the gain for frequency band k depends on NSR^. (n) which in turn 

depends on the noise power, P^, (n) , and noisy signal power, P/ (n) of the same 
frequency band. We have discovered an improvement on this concept whereby 
G, (n) is computed as a function noise power and noisy signal power values from 
multiple frequency bands. According to this improvement, G^. (n) may be computed 
using one of the following methods: 



(1): 



ll-W,{n)^M,NSR, in) 



// = 0.7.27.... 



G,(/zj 



G^ (/I - 1) 



(l-l) 



n = 12. 



T-LT+l 



2T-L... 



// = o,r.2r.... 



Gi(n) = il-VK(/j)^ 



(1.2 



G,(/!-l) 



It = 1.? .... 



r-i.r + i 21-1.... 
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,.w,in)^^ ' . = o,r,2r.... 



Sm,P/(«) (1.3) 
. /i = l,2,...T-l,r + l 27-1,... 



E^*^v(") , /J = 0,7'.27,. 



, n = i,2,...,r- 1,7 + 1 27-1,... 



G,0i-1) 

Our preferred embodiment uses equation (1.4) with M, determined usins the same 
table given above. 

Methods described by equations (1.1)-(1.4) all provide smoothing of the input 
signal spectrum and reduction in variance of the gain factors across the frequency 
bands. Each method has its own particular advantages and trade-offs. The first 
method (1,1) is simply an alternative to smoothing the gains directly. 

The method of ( 1.2) provides smoothing across the noise spectrum only while 
(1.3) provides smoothing across the noisy signal spectrum only. Each method has its 
advantages where the average spectral shape of the corresponding signals are 
maintained. By performing the averaging in (1.2), sudden bursts of noise happening in 
a panicular band for very short periods would not adversely affect the estimate of the 
noise spectrum. Similarly in method (1.3), the broad spectral shape of the speech 
spectrum which is generally smooth in nature will not become too jagged in the noisy 
signal power estimates due to, for instance, changing pitch of the speaker. The 
method of (1.4) combines the advantages of both (1.2) and (1.3). 

There is a subtle difference between (1.4) and (1.1). In ( 1.4). rhe averaging is 
performed prior to determining the NSR rano. In (1.1 ). the NSR values are computed 
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first and then averaged. Method (L4) is computationally more expensive than (LI) 
but performs better than (1.1). 
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Those skilled in the art will recognize that preceding detailed description discloses 
the preferred embodiments and that those embodiments may be altered and modified 
without departing from the true spirit and scope of the invention as defined by the 
accompanying claims. For example, the numerators and denominators of the ratios 
shown in this specification could be reversed and the shape of the curves shown in 
Figures 5, 7 and S could be reversed by making other suitable changes in the 
algorithms. In addition, the function blocks shown in Figure 3 could be implemented 
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in whole or in part by application specific integrated circuits or other forms of logic 
circuits capable of performing logical and arithmetic operations. 
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What is claimed is: 

1 . In a communication system for processing a conmiunication signal derived 
from speech and noise, apparatus for enhancing the quality of the communication signal 
comprising: 

means for dividing said communication signal into a plurality of frequency 
band signals; and 

a calculator generating a plurality of power band signals each having a 
power band value and corresponding lo one of said frequency band signals, each of said 
power band values being based on estimating over a time period the power of one of said 
frequency band signals, said time period being different for at least two of said frequency 
band signals, calculating weighting factors based at least in part on said power band 
values, altering the frequency band signals in response to said weighting factors to 
generate weighted frequency band signals and combining the weighted frequency band 
signals to generate a communication signal with enhanced quality. 

2. Apparatus, as claimed in claim L wherein said calculator comprises a 
memory storing variables having values related to said rime periods which are different 
for at least two of said frequency band signals and wherein said calculator uses said 
variables during said estimating. 

3. Apparatus, as claimed in claim 2, wherein said calculator detects voice 
activity by generating a first signal indicating the probability that said communication 
signal is derived at least in pan from speech and wherein said calculator is responsive to 
said first signa: .WiO herein the values of said \ ariaoies vary depending on the vaiu-j oi" 
said first signal. 
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4. Apparatus, as claimed in claim 3, wherein said power band signals comprise 
noise power band signals each having a noise power band value for one of said 
frequency band signals, each of said noise power band values being based on esdmating 
over a time period the power of noise in one of said frequency band signals, said time 
period being different for at least two of said frequency band signals, wherein said first 
signal has a first value indicating a first probabihty that said communication signal is 
derived at least in part from speech, a second value indicating a second probability 
greater than said first probability that said communication signal is derived at least in 
part from speech and a third value indicating a third probability greater than said second 
probability that said communication signal is derived at least in part from speech, and 
wherein said noise power band values remain substantially constant at least when said 
first signal has said third value. 

5. Apparatus, as claimed in claim 1, and wherein said calculator generates a 
dropout signal in the event that at least one characteristic of said communication signal 
has a defined attribute and wherein said calculator changes the rate at which said power 
band values are allowed to change during the presence of said dropout signal. 

6. Apparatus, as claimed in claim 5, wherein said calculator terminates said 
dropout signal after a predetermined time period. 

7. Apparatus, as claimed in claim 6, wherein said one characteristic comprises 
power of at least one of said frequency band signals. 

8. Apparatus, as claimed in claim 5, wherein said calculator generates a new 
environment signal in the event that said communication signal is detected at the 
beginning of a caii or in the event that said dropout signai has been terminated and 



wo 01/73760 



PCT/USOl/06889 



wherein said calculator changes the rate at which said power band values are allowed to 
change during the presence of said new environment signal. 

9. Apparatus, as claimed in claim 8, wherein said calculator terminates said new 
environment signal after a predetermined time period, 
5 10. Apparatus, as claimed in claim 1, wherein said means for dividing fomis a 

portion of said calculator. 

11. Apparatus, as claimed in claim 1, wherein said calculator comprises a digital 
signal processor 

12. Apparatus, as claimed in claim L wherein said calculator generates a new 
10 environment signal in the event that said communication signal is detected at the 

beginning of a call or in response to at least one characteristic of said communication 
signal having a defined attribute and wherein said calculator changes the rate at which 
said power band values are allowed to change during the presence of said new 
environment signal 

15 13. Apparatus, as claimed in claim 12. wherein said calculator terminates said 

new environment signal after a predetermined time period. 

14. Apparatus, as claimed in claim 3. wherein said communication signal 
defines a variable pitch due to said speech, wherein said system further comprises a pitch 
period detector, wherein said calculator generates a new environment signal in the event 

20 that said pitch period is unsteady and the value of said first signal is greater than a 

predetermined minimum, and wherein said calculator changes the rate at which said 
power band values are allowea to change during the presence of said new environment 
signal. 
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15. In a communication system for processing a communicaiion signal derived 
from speech and noise, a method of enhancing the quality of the communicaiion signal 
comprising: 

dividing said communication signal into a plurality of frequency band 

signals; 

generating a plurality of power band signals each having a power band 
value and corresponding to one of said frequency band signals, each of said power band 
values being based on estimating over a time period the power of one of said frequency 
band signals, said time period being different for at least two of said frequency band 
signals; 

calculating weighting factors based at least in part on said power band 

values; 

altering the frequency band signals in response to said weighting factors to 
generate weighted frequency band signals: and 

combining the weighted frequency band signals to generate a 
conrununication signal with enhanced quality. 

16. A method, as claimed in claim 15, and funher comprising storing variables 
having values related to said time periods which are different for at least two of said 
frequency band signals and using said variables during said estimating. 

17. A method, as claimed in claim 16, and further comprising generating a first 
signal indicating that said communication signal is derived at least in part from speech 
and wherein the values of said variables vary depending on the value of said first signal. 

IS. A methoj. cLimcG in claim i"^. wherein said power band signals comprise 
noise power band signals each having a noise power band value for one of said 
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frequency band signals, each of said noise power band values being based on estimating 
over a time period the power of noise in one of said frequency band signals, said time 
period being different for at least two of said frequency band signals, wherein said first 
signal has a first value indicating a first probability that said communication signal is 
5 derived at least in part from speech, a second value indicating a second probability 

greater than said first probability that said communication signal is derived at least in 
part from speech and a third value indicating a third probability greater than said second 
probability that said conmiunication signal is derived at least in part from speech, and 
wherein said noise power band values remain substantially constant at least when said 
10 first signal has said third value. 

19. A method, as claimed in claim 15, and further comprising: 

generating a dropout signal in the event that at least one characteristic of 
said communication signal has a defined attribute; and 

changing the rate . at which said power band values are allowed to change 
15 during the presence of said dropout signal. 

20. A method, as claimed in claim 19, and further comprising terminating said 
dropout signal after a predetermined time period. 

21. A method, as claimed in claim 20, wherein said one characteristic comprises 
power of at least one of said frequency band signals. 

20 22. A method, as claimed in claim 19, and further comprising: 

generating a new environment signal in the event that said communication 
signal is detected at the beginning of a call or in the event that said dropout signal has 
been lenTiinaied: and 
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changing the rate at which said power band values are allowed to change 
during the presence of said new environment signal. 

23. A method, as claimed in claim 22, and further comprising terminating said 
new environment signal after a predetermined time period. 

24. A method, as claimed in claim 15, and further comprising: 

generating a new environment signal in the event that said communication 
signal is detected at the beginning of a call or in response to at least one characteristic of 
said communication signal having a defined attribute; and 

changing the rate at which said power band values are allowed to change 
during the presence of said new environment signal. 

25. A method, as claimed in claim 24, and further comprising terminating said 
new environment signal after a predetermined time period. 

26. A method, as claimed in claim 17, wherein said communication signal 
defines a variable pitch due to said speech and wherein said method further comprises: 

detecting the period of said pitch: 

generating a new environment signal in the event that said period of said 
pitch is unsteady and the value of said first signal is greater than a predetermined 
minimum; and 

changing the rate at which said power band values are allowed to change 
during the presence of said new environment signal. 

27. In a communication system for processing a communication signal derived 
from speech and noise, apparatus for enhancing the quality of the communication signal 
compnsing: 
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means for dividing said communication signal into a plurality of frequency 
band signals; and 

a calculator generating a plurality of power band signals each having a 
power band value and corresponding to one of said frequency band signals, generating a 
dropout signal in the event that at least one characteristic of said communication signal 
has a defined attribute, changing the rate at which said power band values are allowed to 
change during the presence of said dropout signal, calculating weighting factors based at 
least in part on said power band values, altering the frequency band signals in response 
to said weighting factors to generate weighted frequency band signals and combining the 
weighted frequency band signals to generate a communication signal with enhanced 
quality. 

28. In a communication system for processing a conmiunication signal derived 
from speech and noise, a method of enhancing the quality of the communication signal 
comprising: 

dividing said communication signal into a plurality of frequency band 

signals; 

generating a plurality of power band signals each having a power band value 
and corresponding to one of said frequency band signals; 

generating a dropout signal in the event that at least one characteristic of said 
communication signal has a defined attribute; 

changing the rate at which said power band values are allowed to 
change during the presence of said dropout signal: 

caicuiaimg vveigiiting factors based at least in pan on said power band 

values; 
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altering the frequency band signals in response to said weighting 
factors to generate weighted frequency band signals; and 

combining the weighted frequency band signals to generate a 
communication signal with enhanced quality. 

29. In a communication system for processing a communication signal derived 
from speech and noise, apparatus for enhancing the quality of the communication signal 
comprising: 

means for dividing said communication signal into a plurality of frequency 
band signals; and 

a calculator generating a plurality of power band signals each having a 
power band value and corresponding to one of said frequency band signals, generating a 
new environment signal in the event that said conmiunication signal is detected at the 
beginning of a call or in response to at least one characteristic of said communication 
signal having a defined attribute, changing the rate at which said power band values are 
allowed to change during the presence of said new environment signal, calculating 
weighting factors based at least in pan on said power band values, altering the frequency 
band signals in response to said weighting factors to generate weighted frequency band 
signals and combining the weighted frequency band signals to generate a 
communication signal with enhanced quality. 

30. In a communication system for processing a communication signal denved 
from speech and noise, a method of enhancing the quality of the communication signal 
comprising: 

dividing said communication signal into a plurality ot frequencx bund 

signals: 
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generating a plurality of power band signals each having a power band value 
and corresponding to one of said frequency band signals; 

generating a new environment signal in the event that said communication 
signal is detected at the beginning of a call or in response to at least one characteristic 
of said communication signal having a defined attribute; 

changing the rate at which said power band values are allowed to 
change during the presence of said new environment signal; 

calculating weighting factors based at least in part on said power band 

values; 

altering the frequency band signals in response to said weighting 
factors to generate weighted frequency band signals: and 

combining the weighted frequency band signals to generate a 
communication signal with enhanced quality. 
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FIG. 6 
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An improved noise suppression system (400) which performs speech quality enhancement upon speech-plus-noise 
signal available at the input (205) to generate a clean speech signal at the output (265) by spectral gain modification. The 
noise suppression system includes a background noise estimator (420) which generates and stores an estimate of the back- 
ground noise power spectral density based upon pre-processed speech (215), as determined by the detected minima of the 
post-processed speech energy level. This post-processed speech (255) may be obtained directly from the output of the noise 
supression system, or may be simulated by multiplying the pre-processed speech energy (225) by the channel gain values of 
the modification signal (245). The channel gain controller (240) produces these individual channel gain values for applica- 
tion to both the channel gain modifier (250) and the background noise estimator (420). Each individual channel gain value 
is selected as a function of (a) the channel number, (b) the current channel SNR estimate, and (c) the overall average back- 
ground noise level. The technique of implementing post-processed signal to generate the background noise estimate (325) 
provides a more accurate measurement of the background noise energy since it is based upon much cleaner speech signal. 
As a result, the present invention performs acoustic noise suppression in high ambient noise backgrounds with significant- 
ly less voice quality degradation. 



